Skyrmions and Anomalous Hall Effect in a Dzyloshinskii-Moriya Spiral Magnet 
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Monte Carlo simulation study of a classical spin model with Dzylosinskii-Moriya interaction and 
the spin anisotropy under the magnetic field is presented. We found a rich phase diagram containing 
the multiple spin spiral (or skyrme crystal) phases of square, rectangular, and hexagonal symmetries 
in addition to the spiral spin state. The Hall conductivity Oxy is calculated within the sd model 
for each of the phases. While o^y is zero in the absence of external magnetic field, applying a field 
strength H larger than a threshold value He leads to the simultaneous onset of nonzero chirality and 
Hall conductivity. We find He = for the multiple spin spiral states, but He > for a single spin 
spiral state regardless of the field orientation. Relevance of the present results to MnSi is discussed. 

PACS numbers: 75.10.Hk, 75.10.Jm 



Spiral magnets due to the Dzyloshinskii-Moriya (DM) 
interaction are attracting recent intense attention. MnSi 
is one such prominent example of a metallic chiral ferro- 
magnet whose magnetic [1-3] and transport [4-6] proper- 
ties have been thoroughly studied over several decades. 
Also the spiral magnet is a promising laboratory to study 
the important issue of the relationship between the non- 
trivial spin textures and the quantum electronic trans- 
port, since the non-coUinear spin configurations can be 
easily realized in such magnets. Indeed, Lee et al. found 
the novel anomalous Hall effect (AHE) in a certain re- 
gion of the pressure-temperature phase diagram of MnSi 
recently [5, 6]. Several calculations have pointed out a 
close connection between the spin chirality derived from 
the non-coUinear spin ordering, the Berry phase associ- 
ated with it, and the anomalous Hall effect which arises 
as a consequence[7-10]. 

The purpose of this paper is twofold. First, we ad- 
dress the possible phase diagram of a spiral magnet from 
a microscopic model that includes the DM interaction 
and spin anisotropy. The phase diagram in the plane of 
spin anisotropy-magnetic field including the multiple spin 
spiral or Skyrmion crystal phase of square, rectangular, 
and hexagonal symmetries is revealed. The anomalous 
Hall conductivity is calculated for a model of conduction 
electrons coupled to the spins via the double exchange 
interaction, and the relevance of Skyrmion spin texture 
on the Hall transport is studied. In this way, a clear and 
consistent connection is drawn between the underlying 
topological spin structure, and its manifestation in an 
anomalous Hall transport. 

A continuum Hamiltonian written by Bak and Jensen 
for a prototypical chiral magnet MnSi some years ago[ll] 
is adapted to a lattice spin model consisting of the ferro- 
magnetic exchange (J), DM interaction (K), anisotropy 



{Al and A2), and the Zeeman energy (H): 

Hs = — J Sr • (Sr-ti + Sr-fiJ + Sr+z) 
r 

-if ^ (Sr X SrJ^ -X + SrX S^+y -y+SyX Si^£ • z) 
r 

r 

-A,Y,[s^s^^+sysl^+s^^s^,) - H . ^ Sr. 

r r 

(1) 

We take a cubic lattice structure rather than the full B20 
lattice structure of MnSi[12]. 

Without the anisotropy terms and the magnetic field, 
Hs favors the spiral spin (SS) ground state in which the 
spiral spins lie in a plane orthogonal to the propagation 
vector k: Sr - Ske'"""" + S^e^^'^^ Sk-k = 0. The k 
stretches along the [111] ([11] in 2D) direction with the 
magnitude k = |k| fixed by tanfc = [K / J){\ / ^/d) , d be- 
ing the spatial dimension. The pitch of the helix can be 
very long compared to lattice spacing due to the small 
K/ J ratio in a material. Under high pressure there is also 
experimental evidence for the realization of a phase with 
multiple ordering vectors [3]. Throughout this paper we 
denote such a multiple-k spin structure as Skyrme crys- 
tal (SC). The reason for the nomenclature will become 
transparent shortly. 

Monte Carlo (MC) simulated annealing procedure was 
employed to work out the ground states for varying 
anisotropy strengths (yli,^2) and the field strength H 
oriented along the z-direction: H = Hz. Some simplifi- 
cations were made to save the computational cost. First, 
a 2D rather than the 3D lattice was used. Because the 
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realistic modulation period is very large and difficult to 
simulate, we also choose the ratio K/J = 27r/6 (Here- 
after we will take J = 1) which would give k = 27r/6 
in 2D without the anisotropy. Calculations were mostly 
carried out for 18 x 18 lattice, with occasional checks on 
a 30 X 30 lattice to ensure consistency. 2 x 10^ MC steps 
were used at each temperature in the annealing process. 
It turns out that the same ground state is found over a 
widely different choices of ^i, and here we present all 
the results for Ai = 0.5 without loss of generality. Once 
the ground state has been obtained for a given A2 and 
H, we analyze its structure by making the Fourier trans- 
form (Sk) = X]r(^r)e~**'''" of the averaged MC configu- 
rations (Sr), and looking at the intensity profile |(Sk)p. 
A sharper spectral feature is obtained in this way than by 
taking the average of the individual intensities, (|SkP). 

A pair of sharp peaks are obtained at ±k for a single 
SS, where k = (fc, k). Each spin spiral has a right-handed 
helicity consistent with the sign oi K > 0. When A2 is 
sufficiently large, one finds the SC phase given as the su- 
perposition of two pairs of spin spirals, with modulation 
vectors at (±A;,0) and (0, ±fc). Depending on whether 
the Bragg intensities |(Sk)P are the same or different for 
the two pairs, we denote them as SCi (non-identical) or 
SC2 (identical). In practice SCi is fragile, occupying only 
a tiny fraction of the phase diagram. A third class of SC 
state is found when the field strength in exceed of a cer- 
tain threshold value He is applied to the SS phase. This 
state, denoted SC/i, is characterized by three sets of mod- 
ulation vectors which are related by 120° rotations. We 
later point out that this phase, and SCi and SC2 phases 
under nonzero magnetic field, also carries a nonzero uni- 
form spin chirality and anomalous Hall conductivity. The 
low temperature schematic phase diagram as obtained 
from examination of the Bragg intensity |(Sk)P is shown 
in Fig. 1. 

Although the details of the phase transitions among 
the phases deserves further, more careful study in the 
future work, we can make some qualitative statements 
based on symmetry. The transition from SS to SCh and 
SS to SCi are necessarily first-order, while SCi ^ SC2 
can be second-order. In going from SCh or SC2 to the 
fully spin-polarized (SP) state at high field, the intensity 
of Bragg peaks gradually diminishes, leaving only the 
k = peak in the SP phase. 

The projection of spin patterns onto the xy plane is 
displayed in Fig. 2 for the three Skyrme crystal states. 
We introduce the local chirality Xr at the lattice site r as 

STTXr = Sr • (814^ X Si4^) + Sr • (S,^ X S,^). (2) 

It is well known that a single localized Skyrmion would 
give the uniform chirality x = X^r Xr equal to unity. 
Plots of Xr in Fig- 2 clearly display the presence of 
Skyrmion (bright) and anti-Skyrmion (dark) regions. 
The skyrmion density map is largely one-dimensional for 



SCi, checkerboard- like for SC2, and hexagonal for SCh- 
The uniform chirality x was zero for the states obtained 
at zero field. 

Expanding Sr = X^k Ske* '', the uniform chirality 

reads x = J2k^M ^ki+ks ' x Sk^ cos(A;i:r + ksy). Our 
calculation explicitly shows that this quantity is zero 
without the magnetic field. With the field on there ap- 
pears the k = component of spin. So oc H. One 
might think that a non-zero chirality results from x ~ 
H • Sk X S_k, but since Sk x S_k ^ ik and thus 
Sk X S_k + S_k X Sk = 0, the chirality contribution from 
a given spin spiral is zero. Rather, we find that the uni- 
form chirality is induced in the SCi and SC2 phases due 
to the appearance of additional Bragg peaks at S±k^±k„j 
ka; = (A;,0) and k^ = (0, fc), and the nonzero value 
of the triple product such as Sk^.^^^ • Sk^ x Sk„ . Why 
the Sk^+k„ component appears with H can be under- 
stood from a Ginzburg-Lanau theory which should con- 
tain, on symmetry grounds, a quartic coupling such as 
(Ski-Sk2)(Sk3-Sk4). With the appearance of Sq ^ H, the 
coupling becomes (Sk^ •Sk„)(Sk^_^(.^ -H). Since Sk^-Sk^ 
is already nonzero in the SC phase, it implies the ap- 
pearance of a linear term in S'k^^jc^ which will lead to 
the condensation of S'k^-fjc order and x proportional to 
H. The numerical results shown in Fig. 3 confirm the 
linear growth of x with H. For SC/i, the imiform chi- 
rality is found to arise from the nonzero triplet product 
X ^ Ski ■ Ska X Ska ' where the three independent modu- 
lation vectors form ki -|- k2 -|- ka = 0. 

For field oriented away from the 2;-axis, the Ginzburg- 
Landau argument would continue to predict the presence 
of a linear- dependence of the chirality in the SCi and 
SC2 phases and lack therefore for the SS phase. Both 
these predictions are confirmed by the numerical calcu- 
lations. Some check was made on the finite-temperature 
behavior, using the thermal-averaged intensity |(Sk)P 
was used as a measure of magnetic ordering. The largest 
intensity in |(Sk)P among all vahies of k decreased con- 
tinuously until it reaches zero at the critical temperature. 
There was no clear evidence for the existence of an inter- 
mediate phase. 

Having studied the magnetic phase diagram, we turn 
to the coupling of the local moments to the conduction 
electrons that would result in the anomalous Hall effect. 
We adopt the sd Hamiltonian 



H = -t^ C+^CrV - A ^ Sr • C+„(cr)a/3Cr/3, (3) 
rr'cT r 

with the moment distribution {Sr} obtained from previ- 
ous MC calculation. The intrinsic anomalous Hall con- 
ductivity axy is calculated from 
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(4) 

expressed in units of e'^/h. The sum J2m^n extends over 
all non-identical pairs of singlc-particlo oigonstatcs of Eq. 
(3), fm is the Fermi function 1/(6^^"* + 1), and gives 
the number of lattice sites. and Jy are the current 
operators. The same temperature T is used both in the 
Monte Carlo generation of sample spin configurations, 
and in evaluating axy Thermal average {(Txy) was taken 
over 100 MC-generated spin configurations. The relax- 
ation rate rj = 0.1, and the sd coupling A = 1 were used 
in the calculation, with similar results at other parameter 
choices. Since the relativistic spin-orbit coupling term is 
absent in Eq. (3), the nonzero Hall conductivity we will 
report below can only be due to the topological Berry 
phase effects. 

We find {axy) = for all the spin configurations at all 
temperatures when H = 0. For finite fields, the onset of 
nonzero {(Txy) and nonzero x coincided almost perfectly, 
confirming the earlier theoretical anticipation that an un- 
conventional anomalous Hall conductivity arises as a con- 
sequence of uniform chirality x being nonzero[7, 8, 10]. 
We also verified that the relation between {axy) and x are 
linear when both quantities are sufficiently small. Such 
a close tie between x and {(Txy) suggests that measure- 
ment of axy in a given material can be used as an effective 
probe of the underlying spin structure. For instance, lin- 
ear rise of axy with the magnetic field would be consistent 
with the Skyrme crystal spin structure, but not with the 
spiral spin structure. 

The field dependence of uniform magnetization M, 
uniform chirality X) and {axy) are displayed in Fig. 3 for 
several temperatures. Under the full polarization M = 1, 
{axy) and x naturally must go to zero, resulting in the 
characteristic dome shape of the {axy) curve. A roughly 
linear relationship axy ocH ocM exists in both SCi and 
SC2 for small M. The slope axy/H decreases gradually 
as the temperature rises. The low-temperature axy can 
reach up to « lOe^/h, similar to the value reported in 
the model calculation by Metalidis and Bnmo[9]. 

Now we discuss the possible relevance of our results 
to MnSi. Several mechanisms for de-stabilizing the SS 
order in a chiral magnet have been proposed[13-16]. 
Tewari et al. argued for the existence of a chiral liq- 
uid phase on the high pressure side of MnSi, i.e. a 
non- magnetic phase with nonzero (^r). V^r being the chi- 
ral order parameter V'r = Sr • V x Sr. Provided that 
one can associate their pressure variable with A2 in our 
model, we find the skyrme crystal phase characterized by 
Xr = Sr- (i9a:Sr X OySr) in place of their chiral liquid order. 
Robler et a/. [15] considered a Ginzburg-Landau model 
which supports a Skyrme crystal phase at the intermedi- 
ate temperature range between the low-temperature SS 



and the high-temperature paramagnet. Our model in 
contrast exhibits the skyrme crystal phase at low temper- 
ature, driven by the large anisotropy term comparable to 
J and/or K. 

Several features of the experiments on MnSi might be 
understood within the framework of our model. The first- 
order change with A2 in going from SS to SC phase and 
the re-orientation of the primary ordering vector from 
[11] to [10] direction in our model are consistent with the 
first-order collapse of Tc at the critical pressure and the 
re-orientation of the main Bragg peak from [111] to [110] 
in MnSi, respectively[3] . The appearance of an anoma- 
lous Hall signal on the high pressure side [6] and over a fi- 
nite temperature and field segments of the phase diagram 
might indicate the realization of a topologically nontrivial 
spin structure, according to the close tie between anoma- 
lous Hall conductivity and nonzero spin chirality found in 
our study. It is also tempting to associate the recently ob- 
served hexagonal Bragg spots in the so-called A-phase of 
MnSi[18] with the SC;j found in our spin model. Further 
measurement of the anomalous Hall conductivity in the 
same phase will help clarify the connection. The highly 
nonlinear threshold behavior of axy in the SS phase has 
not yet been clearly resolved in existing measurements of 
ordinary Hall effect by Lorentz force or the anomalous 
Hall effect on MnSi[5, 6]. 

We close with a caveat that the correspondence be- 
tween our theory and experiments on MnSi rests on the 
physical assumption that increasing pressure does lead 
to the increase of A2 over other energy scales. Whether 
the Mn orbitals responsible for the localized moments 
in MnSi will give rise to such an interaction, and if the 
pressure does have an effect on A2 remains a subject for 
future consideration. Nevertheless, the predicted field 
dependence of axy on the underlying spin structure is 
independent of the microscopic spin model being used. 
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FIG. 1: (color online) Low temperature (T = 0.01) phase di- 
agram of the spin model in Eq. (1) with K = VG, Ai = 0.5, 
and H = Hz. Phase boundaries are drawn on the basis of 
MC simulations at a largo number of {A2, H) locations. Spin 
configurations are abbreviated as SS (spiral spin), and SC 
(spin crystal). Full spin polarization (SP) results at high field. 
The spin crystal phases are further classified as SCi (unequal 
Bragg intensities), SC2 (equal Bragg intensities), and SC^ 
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FIG. 2; (color online) A plot of the spin configuration pro- 
jected on the xy plane {Sf,S^) in the three spin crys- 
tal ground states: (a) SCi at {Ai,A2,H) = (0.5,2.0,0.0), 
(b) SC2 at {Ai,A2,H) = (0.5,3.0,0.0), and (c) SCh at 
{Ai,A2,H) = (0.5, 0.0, 2.0). At the bottom left of each figure 
are the plots of the Bragg intensity |{Sk)P showing two (SCi, 
SC2) and three (SC^) sets of modulation vectors. Shown at 
the bottom right are the plots of the local chirality Xr. Bright 
(dark) regions correspond to Skyrmions (anti-Skyrmions). 
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FIG. 3: (color online) Upper panels: Plot of uniform mag- 
netization A'l (blue circle) and uniform chirality x (red tri- 
angle) with varying field strength H , H = Hz, for (a) SS at 
{Ai,A2) = (0.5,0.0), (c) SCi at (^1,^12) = (0.5,2.0), and (e) 
SC2 at {Ai,A2) = (0.5,3.0), at three temperatures T = 0.1 
(top six figures), 0.5 (middle), and 1.1 (bottom). Lower pan- 
els: {(Jxy) averaged over 100 MC-generated spin configura- 
tions are shown in (b), (d), and (f) at the corresponding 



